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S2 '. Abstract 

B 

I, ' We report thermodynamic measurements in a magnetic-field-driven quantum critical 

^ ■ point of a heavy fermion metal, YbRh2Si2. The data provide evidence for an energy scale 

O '. in the equilibrium excitation spectrum, that is in addition to the one expected from the 

slow fluctuations of the order parameter. Both energy scales approach zero as the quantum 
critical point is reached, thereby providing evidence for a new class of quantum criticality. 



Quantum criticality encodes the strong fluctuations of matter undergoing a second order 
phase transition at zero temperature. It underlies the unusual properties observed in a host 
of quantum materials. A basic question that remains unsettled concerns its proper theoreti- 
cal description, which is challenging because the fluctuations are both collective and quantum 
mechanical. One class of theory, based on the traditional formulation of classical critical phe- 
nomena [[T]|, considers the fluctuations of a classical variable - Laudau's order parameter - in 
both spatial and temporal dimensions [|2l[3l|4l[51. The slowing down of the order-parameter fluc- 
tuations accompanies the divergence of a spatial correlation length; at each value of the tuning 
parameter, the equilibrium many-body spectrum contains a single excitation energy scale, which 
vanishes at the quantum critical point (QCP) fS]. An unconventional class of theory lEllHlIll, by 
contrast, is inherently quantum mechanical; it explicitly invokes quantum entanglement effects, 
which are manifested through vanishing energy scale(s) that are in addition to the one associ- 
ated with the slowing down of order-parameter fluctuations. The nature of quantum criticality 
can therefore be experimentally elucidated by determining whether a single or multiple energy 
scale(s) vanish as the QCP is reached. 

We consider the heavy-fermion metal YbRh2Si2 (YRS), and show that multiple energy 
scales vanish as its QCP is approached and, in addition, suggest that critical electronic modes 
co-exist with the slow fluctuations of the magnetic order parameter. A direct way to probe the 
intrinsic energy scales in the equilibrium spectrum near a QCP is to measure thermodynamic 
properties. Another approach is to measure the fluctuation spectrum in equilibrium, for exam- 
ple by inelastic neutron scattering experiments. Such equilibrium methods are in contrast to 
transport experiments, which are influenced by electronic relaxational properties, especially for 
anisotropic and multi-band systems. 

As extraction of critical energy scales requires measurements through fine steps of the con- 
trol parameter, which is nearly impossible for inelastic neutron scattering, we report here mea- 



surements of thermodynamic properties of YRS across its magnetic QCP. 

We choose to work with the tetragonal heavy fermion compound YRS as it presents a clean 
and stoichiometric material that is well characterized ifTOl . In the absence of an external mag- 
netic field, YRS shows very weak antiferromagnetic (AF) order at T/v = 70 mK with an ordered 
moment of only ~10~^/iB/Yb [fTT|. A small magnetic field {H±c ~ 0.06 T, for the field applied 
within the easy ab plane, and ifyc = 0.66 T, along the hard c axis) suppresses the transition 
temperature and accesses the QCP [[T2|. The ability of using such a small magnetic field to 
access the QCP makes YRS suited for our purpose; the determination of energy scales requires 
scanning across the phase transition, and an external magnetic field can be tuned with relative 
ease and continuously. Hall effect measurements lfT3l on YRS have shown a large and rapid 
crossover in the Hall constant at a temperature-dependent magnetic field away from the anti- 
ferromagnetic transition. In the zero-temperature limit, this crossover extrapolates to a jump 
across the QCP, which has been interpreted as a large change of the Fermi surface volume. This 
represents yet another advantage of measuring the thermodynamic properties in YRS, as they 
can be compared with their transport counterparts. 

We measured the isothermal linear magnetostriction 9 In L/dH, where L is the length along 
the [110] direction within the tetragonal ab plane and the magnetic field H is applied along the 
same direction (H ±c). Fig. 1 shows the magnetostriction as a function of the magnetic field, at 
temperatures ranging from 0.02 K to 0.8 K. For temperatures below 0.075 K, a clear disconti- 
nuity is observed when suppressing the AF order by a critical magnetic field. At T > 0.075 K, 
it is seen that, for small magnetic field, the isothermal magnetostriction linearly depends on the 
magnetic field, as is the case in typical metals [[T4|. Beyond a crossover field, however, there 
is a change to a high field region with a different slope. The crossover field decreases as the 
temperature is reduced. 

To understand this crossover, we compare it with the field-dependent isothermal behavior of 



other thermodynamic and transport quantities. Fig 2A illustrates the similarity of the crossover 
in the magnetostriction with that seen in the field-dependent isothermal Hall resistivity pn (mea- 
sured with H\\c). The Hall coefficient was described [[T3l by an empirical crossover function of 
the form /(if, T) = y42-(^2-^i)/[l + (^/^o)1; thecrossoverfieldscalei7o(r) is equivalent 
to an energy scale T*{H). We have analyzed the magnetostriction data, as well as the existing 
magnetization data (H ± c) [fTSl [T6l . with the same crossover function. Note that no corre- 
sponding anomalies can be resolved in the magnetization data for if || c lfT6l . which is almost 
linear in H. The solid curves in Fig. 1 and Fig. 2A correspond to fits of A[iio], M = M + xH 
and the Hall resistivity pn- Fig. 2B shows the three sets of Hq{T), obtained from such fits. 
Their overlap represents a key conclusion of the present work; it suggests that they define one 
energy scale T*(H). This scale is seen to be distinct from either the transition temperature (T/v) 
for the magnetic ordering at H < He or the scale {Tlfl) for the establishment of the Landau 
Fermi liquid state at if > He. For all three quantities, the width of the crossover extrapolates to 
zero at T = 0, implying that the differentials of the magnetostriction, magnetization, and Hall 
resistivity have a jump in the zero-temperature limit (supporting online text). 

The results raise the important question of the causal relationship between the thermody- 
namic and electronic transport properties. One might argue ifTTl that the Hall-effect evolution 
as a function of the magnetic field lfT3l is caused by the Zeeman splitting of the Fermi surface 
induced by the magnetization (and reflected in the magnetostriction). However, the magnetiza- 
tion only displays a smeared kink, and the corresponding Fermi surface change would at most 
produce a smeared kink in the Hall coefficient evolution; such a kink is too weak compared 
to the smeared jump seen experimentally. Moreover, along the c— axis, even such a smeared 
kink feature is absent in the magnetization vs the magnetic field. Instead, it is more natural to 
attribute the nonanalyticities in both the magnetostriction and magnetization as thermodynamic 
manifestations of the large Fermi surface jump caused by /—electron localization. 



To explore this issue further, we have also studied the longitudinal magnetoresistivity. Fig. 3 
shows the electrical resistivity, p, as a function of the magnetic field (if ± c), at various tem- 
peratures. The broadened step-like decrease, observed at all temperatures, corresponds to the 
crossover observed in the other properties. Indeed, as shown in Inset A, the crossover fields 
determined from the minima of the derivative dp/dH (cf. Inset B) fall on the same T*{H) 
line determined from the magnetostriction, magnetization, and Hall effect. In addition. Inset B 
shows that the width of the crossover decreases as temperature is lowered. A detailed analysis 
shows that the crossover width goes to zero in the zero-temperature limit (supporting online 
text), implying a jump in the residual resistivity across the magnetic QCP. This is in accordance 
with the theoretical expectations lfT8l[T9l associated with an /-electron localization transition. 

Fig. 3, Inset A, also shows the temperature scale as a function of field, extracted from the 
peak in the T-dependence of the differential susceptibility Xac = dM/dH; the latter, observed 
earlier [|20ll . necessarily accompanies the smeared kink behavior in the isothermal M vs H. It is 
clearly seen that this scale too falls on the same T*{H) line. 

Our results shed light on the overall phase diagram of this clean stoichiometric quantum crit- 
ical material. NMR measurements f2\\. while signaling the dominance of AF fluctuations in the 
quantum critical regime, have also revealed enhanced ferromagnetic fluctuations: the Korringa 
ratio, S = 1/TiTK'^, is small - of the order of O.lS'o, where Sq is the corresponding ratio for 
non-interacting electrons. Further evidence for enhanced ferromagnetic fluctuations has come 
from magnetization measurements EOl : the Wilson ratio, Rw = Ti:'^k\/{piQpiljj) x Xo/7oj with 
p,eff = 1.4 ubI^^ [l22|. is strongly enhanced for an extended region of the phase diagram; it is 
already large (~ 20) for magnetic fields of a few Teslas, and further increases as the field is re- 
duced towards He. Therefore, it could be tempting to consider the q ~ magnetic fluctuation as 
the dominant critical fluctuation [17], especially since a conventional ferromagnetic QCP would 
yield a Griineisen exponent ll23l of 1/ zu = 2/3, close to what is observed in YRS [24]. This 



picture is problematic for a number of reasons, however. First, neither 3D nor 2D ferromagnetic 
spin fluctuations can generate the fractional exponent observed in the temperature dependence 
of the uniform spin susceptibility ll20l . Second, ferromagnetic spin fluctuations would lead to 
a divergent 1/Ti (~ 1/T^, with x = 1/3, 1/2 for 3D and 2D cases, respectively), that is in 
contrast to the observation of 1/Ti ~ const, when the NMR measurement field is extrapolated 
to the quantum critical regime [fTTI . Third, as ferromagnetic spin fluctuations are inefficient in 
affecting charge transport, this picture contradicts the observation of a nearly if —independent 
ratio A/x^ that accompanies a strongly if —dependent A and x ll20l . Here, A is the coefficient 
of the T^ component of the resistivity. 

The data presented here show that the uniform (i.e., q = 0) magnetization depends on the 
same underlying physics as that for the charge transport. Since the transport is dominated 
by large q fluctuations, the results imply that the q = magnetic fluctuations are a part of 
overall fluctuations in an extended range of wavevector scales. It is then more natural to assume 
that the dynamical spin susceptibility at different wavevectors obeys the same form ll25l [51 
as observed in another prototypical quantum critical heavy fermion metal, CeCus.gAuo.i ifTTl : 
x(q,T,tu) ~ [0q + T"M/(cj/T)]~\ At the QCP, the Weiss field at the antiferromagnetic 
wavevector (q = Q) vanishes: 0q = 0. At the same time, and unlike for CeCus.gAuo.i, 0q=o 
is very small in YRS. Based on the saturation scale seen in the temperature dependence of the 
uniform magnetic susceptibility ll20ll and the NMR Knight shift data ||2T1 near He, we estimate 
0q=o to be of order 0.3 K [|26l . When q moves away from either or Q, 0q increases to 
the order of the RKKY interaction or bare Kondo scale (about 25 K for YRS [fTOlD . This is 
illustrated in Fig. 4. The enhanced uniform magnetic susceptibility, the concomitant enhanced 
Wilson ratio [|20l . as well as the small S = 1/TiTA'^, naturally follow from this picture. 
Moreover, both Xq=o and Xq=Q scale similarly with H and the observation that A/x^ is nearly 
if -independent is in fact a manifestation of an if -independent A/xq- All these lead to the 



conclusion that the origin of the T* line lies in an electronic slowing down and, for YRS, the 
strong q = fluctuations happen to be a consequence of the latter as well. 

We now turn to more detailed theoretical implications of our results. Our measurements es- 
tablish that the energy scale T* is associated with the equilibrium many-body spectrum (which 
alone determines thermodynamics). Moreover, this scale is distinct from the Landau Fermi liq- 
uid scale, Tlfl, as physical quantities manifest rather different behavior across the two scales 
(supporting online text). Finally, both of these scales vanish at the QCP. These findings con- 
tradict the conventional order-parameter fluctuation theory in at least two respects. First, the 
only low-energy scale in that theory is associated with the magnetic slowing down which, for 
H > He, is Tlfl EEmKH. Second, within that theory, a sharp feature in thermodynamics 
and transport quantities might arise near T^ only. 

Our results are instead consistent with magnetic quantum criticality accompanied by the 
destruction of Kondo entanglement. In the form of local quantum criticality |l7l[8]|, a collapse of 
a large Fermi surface as H decreases leads to an added energy scale characterizing an electronic 
slowing down and, in addition, yields a zero-temperature jump in the Hall coefficient and in the 
field-differentials of the thermodynamic quantities. An additional energy scale also exists in the 
"deconfined" quantum criticality scenario for insulating quantum magnets flUl, as well as in its 
extension to itinerant electron systems Il27ll28l that are argued to be relevant to quantum critical 
heavy fermion metals. 



Figure 1: Magnetic field dependence of the magnetostriction of YbRh2Si2. The symbols rep- 
resent the linear coefficient A[iio] = d\nL/dH (where L is the sample length along the [110] 
direction within the tetragonal ab plane) vs H at various temperatures. Note, that A [no] < 
and that the data sets have been shifted by different amounts vertically. The sharp feature in 
the 0.02 K data corresponds to a discontinuity in A (as is more clearly seen in the measured 
length vs H, which shows a change in slope but does not contain any discontinuity), demon- 
strating the continuous nature of the magnetic transition at the critical field of 0.05 T. Similar 
behavior is observed at various different temperatures below 0.075 K, e.g. at 0.03 K, 0.04 K, 
0.05 K and 0.06 K. The solid lines for T > 0.13 K are fits using the integral of the crossover 
function f(H,T) (see text) which reveal a characteristic field scale Hq(T) along which the 
magnetostriction shows a drastic change in slope. 



Figure 2: Energy scales in YbRh2Si2, determined from thermodynamic, magnetic and transport 
measurements. (A) the field dependence {H ± c) of the magnetostriction A[iio], M = M + xH 
(where x = dM/dH), and the Hall resistivity p^ for H \\ c (for the latter, the field values 
have been multiplied by the factor 13.2 which is the anisotropy ratio of the critical field parallel 
and perpendicular to the c-axis), respectively, all at T = 0.5 K; similar behavior is observed at 
other temperatures. For \[iio](H ± c) the sample has po = 0.5 /ificm and He = 0.05 T, for M 
(H ± c) the sample has po = 1-0 ptl^cm and He = 0.06 T, and for ph{H \\ c) the sample has 
Po = 1-0 /ificm and He = 0.66 T. The solid lines correspond to fits using the integral of the 
crossover function f{H,T) (see text). Each data set has been normalized by its initial slope. 
For clarity, the three data sets have been shifted by different amounts vertically, ph — PH,a, 
where pH,a is the anomalous Hall resistivity iTTSl . behaves similarly as pn- We analyzed M{H), 
which is the field derivative of the magnetic free energy contribution (— M x H)\ fitting M{H) 
leads to similar conclusions, though the quality of the fit is somewhat poorer since M vs H is 
not as linear as M at high fields (Fig. S2). (B) the crossover field scale Hq as determined from 
magnetostriction, M and Hall resistivity using the same symbols as in A. It is equivalent to an 
energy scale T*{H). The gray diamonds and triangles represent respectively the Neel ordering 
temperature (T/v) and the crossover temperature (Tlfl) below which the electrical resistivity 
has the Fermi liquid form, p = po + AT'^, as determined from measurements on a single crystal 
with Po = 0.5 pricm and He = 0.05 T. The solid and dotted lines are guides to the eye; for the 
latter, data points outside the plotted field regime have also been used. The horizontal error bars 
represent the fitting error rather than the width of the crossover. 



Figure 3: Longitudinal magnetoresistivity of YbRh2Si2 as p vs H {H 1. c) at various temper- 
atures. The maxima in the 0.03 and 0.07 K data indicate the boundary of the AF ordered state 
(Tn{H = 0) = 0.075 K). The arrows mark the positions of inflection points in p{H). The inset 
A displays the phase diagram, where the gray shaded area represents the range of Hq{T) values 
shown in Fig. 2B. Open yellow triangles mark the positions of the inflection points in the longi- 
tudinal electrical resistivity. The smeared kink behavior in the isothermal M vs H corresponds 
to a peak in the T-dependence of x- The latter has been observed [|20l ; the corresponding peak 
temperature vs H for a sample with po = 0.5 jiVtcm and He = 0.05 T is displayed by the 
dark blue diamonds and found to be consistent with Hq(T). The inset B displays the derivative 
dp/dH vs H at both T = 0.1 K and 0.3 K. Arrows mark the minima, corresponding to the 
inflection points in p{H). 

Figure 4: Sketch of the suggested q dependence of the Weiss temperature 0q, which enters 
the magnetic susceptibility. As H reaches H^ it vanishes at the antiferromagnetic wavevector 
Q, as shown. In addition, 0q has a second minimum at q = 0. Ferromagnetic fluctuations at 
q = remain important as ©q=Q goes to zero. 
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SOM 1. Relationship between T* and Tlfl 

The two scales, T* and Tlfl, are distinct due to the rather different behavior that physical 
properties manifest across these scales. 

Across the T* line, as a function of the magnetic field {H), the behaviors of isothermal 
physical quantities are governed by a zero-temperature jump that is smeared by temperature. 
The physical quantities include d\/dH, x = dM/dH, dpu/dH, and p. The FWHM, defined 
in the main text, for each of these quantities is plotted as a function of temperature in Fig. SI A. 
Each FWHM extrapolates to zero in the T = limit, signaling a jump in the zero-temperature 
limit. For each quantity, the ratio A2/A1 (where A2 and Ai, defined in the main text, are respec- 
tively the low-field and high field values) is plotted as a function of temperature in Fig. SIB. It 
extrapolates to a value significantly below 1 in the T = limit, implying that the extrapolated 
jump at zero temperature is non-zero. The existence of this zero-temperature jump implies that 
none of the four quantities obeys a scaling in terms of T /T* . 

The Tlfl line, on the other hand, represents a very different crossover. This line was already 
extensively discussed, based upon the specific heat and resistivity data, in Ref. 22 of the main 
text. For resistivity, this scale appears in the temperature dependent component. Indeed, dp/dT 
satisfies a scaling in terms of T/Tlfl- 

dp f T \ 

Alternatively, the contrast between the Tlfl line and the T* line can be seen by integrating 
the scaling equation ©, yielding a p{H, T) — p{H, 0) vs. H that always has a smooth crossover 
across the Tlfl line. This is in contrast to the jump associated with the T* line discussed earlier. 

SOM 2. M vs. H 

As shown in Fig. S2, the magnetization itself displays a similar crossover as M = M + xH 
does. Compared to the M case, its linear H dependence in the measured high field regime is 

16 



somewhat less robust, making the fit by the fitting function f{H, T) to be of a slightly lower 
quality. 
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Figure SI: (A) Temperature dependence of the full- width at half-maximum (FWHM) derived 
from fitting transport and thermodynamic measurements across T*{H) by the crossover func- 
tion }\H,T) = A2- {A2 - Ai)/[l + {H/HqY]. For the Hall resistivity pn for H \\ c, the 
values have been divided by 13.2 which is the anisotropy ratio of the critical field parallel and 
perpendicular to the c-axis. (B) Temperature dependence of the ratio A2/A1. 
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Figure S2: Magnetic field dependence of the Magnetization M(H) of YbRh2Si2 at 0.16 K 
(triangles) and 0.5 K (circles). Lines are fits according to the crossover function with Bq = 
0.11 T (0.16 K) and Bq = 0.2 T (0.5 K). 
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